Phonon-mediated negative differential conductance in molecular quantum dots 
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Transport through a single molecular conductor is considered, showing negative differential con- 
ductance behavior associated with phonon-mediated electron tunneling processes. This theoretical 
work is motivated by a recent experiment by Leroy et al. using a carbon nanotube contacted by an 
STM tip [Nature 432, 371 (2004)], where negative differential conductance of the breathing mode 
phonon side peaks could be observed. A peculiarity of this system is that the tunneling couplings 
which inject electrons and those which collect them on the substrate are highly asymmetrical. A 
quantum dot model is used, coupling a single electronic level to a local phonon, forming polaron 
levels. A "half-shuttle" mechanism is also introduced. A quantum kinetic formulation allows to 
derive rate equations. Assuming asymmetric tunneling rates, and in the absence of the half-shuttle 
coupling, negative differential conductance is obtained for a wide range of parameters. A detailed 
explanation of this phenomenon is provided, showing that NDC is maximal for intermediate electron- 
phonon coupling. In addition, in absence of a gate, the "floating" level results in two distinct lengths 
for the current plateaus, related to the capacitive couplings at the two junctions. It is shown that the 
"half-shuttle" mechanism tends to reinforce the negative differential regions, but it cannot trigger 
this behavior on its own. 

I. INTRODUCTION 

The prospect of using molecules as the fundamental building blocks of future nanoelectronics devices is rather 
innovating and exciting from the point of view of potential applications. On the fundamental side, the field of molecular 
electronics also opens new directions for research because of the prominent role of phonon (vibron) excitations in 
electronic transport. The nano-objects which are connected to metallic leads may consists of individual molecules, 
self- assembled monolayers, or conjugated systems such as polymers and carbon nanotubes. The interplay of electron 
transport and molecular vibrations has triggered much interest, and unambiguous signatures of phonons have been 
detected in several experimentsiSiiiiSiSiL^. 

A recent work 7,8 considered electron injection from a scanning tunneling microscope (STM) tip into a carbon 
nanotube. A single-wall carbon nanotube (CNT) is freely suspended over a trench. The STM-tip is located near 
the center of the suspended part of the CNT. A DC-bias voltage V is applied between the substrate and the STM, 
and the current flowing through the STM-tip - CNT - substrate structure is measured at a given tunneling distance, 
controlled by the setpoint current. The motivation for this special geometry is to allow for free internal vibrations to 
occur in the suspended portion of the nanotube, in particular the so-called radial breathing modes (RBM)S. Contrary 
to acoustic modes in such system, this mode has a rather high oscillation frequency, and the authors observed that 
conductance peaks in the current- voltage characteristics of the nanotube were surrounded by phonon side peaks, due 
to emission or absorption of RBM phonons. Indeed, CNT display a sizeable coupling of electrons to RBM modest. 
The contacts to the nanotube being rather resistive, transport in this system is dominated by the Coulomb blockade 
regime, and phonon side peaks were observed around each Coulomb blockade peak in differential conductance plots. 

Interestingly, the authors mention frequent detection of negative differential conductance (NDC) regions. Striking 
NDC features also appear in a very recent work by Sapmaz et al&, in a transport measurement of a suspended CNT, 
where current flows through the CNT, between two contacts at the substrate. In this work, phonon side peaks are 
attributed to longitudinal stretching modes, and the steps in the I(V) characteristics are followed by spikes, thus 
displaying NDC features. The purpose of the present theoretical work is to show that such NDC features can be 
described quite simply using a generic model which consists of a quantum dot with a single orbital level, coupled 
to on-site single phonon mode, and connected to leads by tunnel junctions. Due to the weak tunnel couplings, the 
physics of NDC appear to be a consequence of the transport through small polaron states on the molecule. Notice 
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that polaron formation in CNT is suggested by several worksiS. Similar models have previously been considered in 
the Iiteratureiiii2ii&i!ii5ii£ii!. Here we consider a molecular system - or a nanotube setup - whose tunneling matrix 
elements from the molecule to the leads are asymmetric, which is the case in Ref. Q, and possibly also in Ref. 0. 
We propose that given this asymmetry, NDC can be obtained for a wide range of parameters and it can even be 
quantified by analytical means for strong asymmetry. Furthermore, we argue that our approach for the description of 
NDC is by no means confined to the experimental geometry of Ref. 0, which can be considered as an "experimental 
paradigm" for phonon modulated molecular transport. The present approach should apply to any molecular transport 
in the Coulomb blockade regime where an optical-like phonon mode dominates. Notice that NDC was found in RefUll 
within a symmetric model, using an ansatz for the Franck-Condon factors coupling the leads to the polaron levels. 
Catastrophic NDC was also proposed to occur in Ref. IT3. but due to an additional mechanism. In the framework of 
a single site model, NDC was also found in an adiabatic treatment of the polaron problem^. It has been proposed to 
occur in the case of two competing molecular stated or within a two-site modeli 9 .. Besides electron-phonon coupling, 
other physical mechanisms can lead to NDC-S. 

In molecular electronics, one is tempted to use existing theoretical tools developed in the context of mesoscopic 
physics because the size of the devices allows a coherent description of the transport inside the molecule. When the 
molecule is connected to good contacts, and when the role of Coulomb interaction and electron-phonon interactions is 
reduced, the Landauer formulation of transport combined with a Green's function calculation of the transmission coef- 
ficient serves as a good starting point for computing transport^, but it then neglects the dynamical degrees of freedom 
of the molecule. Phonons can be included in these approaches, analytically 23 , or using numerical calculations 2 ^ 2 ^ 
based on the work of Ref. 12a . Alternatively, the Landauer approach can be generalized to a situation where interac- 
tions are restricted to a finite regioniSiiSiSi using non-equilibrium Green's functions. 

On the other hand, molecular contacts or the contacts to a nanotube are often of poor quality. It is then reasonable 
to think of the molecule as a quantum dot, which is subject to the Coulomb blockade. This point of view was adopted 
by several authorsii^iSii^i^iiSi 2 ^, and is confirmed by experiments. One of the advantages of modecular electronics 
comes from the size of the individual nano-object which is connected to the electrodes. Even at relatively large 
temperatures (a few Kelvins), the dynamics within the molecule/nanotube is fully phase coherent. On the other 
hand, the tunneling rates to the reservoirs are typically small compared to temperature. One can then derive rate 
equations for the electron population of the dot. Here our goal is to use the simplest model which can account for 
the physics of the experiment of Ref. Q- We therefore derive a kinetic equation approach which can account for the 
observed NDC effect, treating the electron-phonon coupling non-perturbatively while still allowing for an intuitive 
understanding of the physics at hand. Furthermore, consequences can be derived on the electron transport as well as 
on the non-equilibrium phonon population. In the present work the local mode is not coupled to any environmental 
degrees of freedom, like phonon modes or electron-hole excitations in the substrate. Although such coupling is probably 
relevant in some experiments, the very high phonon quality factors (Q > 20000) obtained in other worksi justify to 
neglect it as a first step. 

Experimentally, the addition of voltage gates on a molecular transport setup is still challenging but it is possible 8,29 
to approach three metallic probes in a nanometer scale region. Yet, there is also a motivation to study a setup where 
the molecular levels are "floating" instead of being fixed by a gate. Besides the strong asymmetry of the tunneling 
rates, a specific feature of our work is that NDC can occur in such absence of gate voltage, and that the ratio of the 
right/left capacitance plays an important role, as it dictates the location of the molecular levels. 

The structure of the paper is as follows. The model is introduced in Sec. II, and the derivation of the rate 
equations is provided in Sec. III. The general considerations about current transport and phonon occupation numbers 
are discussed in Sec. IV. Numerical results are presented in Sec. V, illustrating the role of the asymmetry of the 
capacitances, of the shift of the dot level, and of the so-called half-shuttle mechanism on the NDC. We conclude in 
Sec. VI. 



II. THEORETICAL MODEL 

Although the potential relevance of our approach applies to the carbon nanotube experiments of Refs. I7l8l we 
refer to the central region between the leads as the quantum dot. This quantum dot is weakly coupled to two 
metallic electrodes by tunnel junctions. In the STM geometry of Refs. lllcSl the electrodes represent the STM-tip 
and the substrate, to which we refer as the 'left' (I) and 'right' (r) electrodes, respectively. Each tunnel junction 
(j = l,r) is characterized by a resistance (Rj) and a capacitance (Cj). While R r and C r are constant for a given 
contacted nanotube sample, Ri and Ci are functions of the tip-tube separation. The capacitances C; and C r are 
sample-dependent. For generality, one can add a gate voltage V g and a gate capacitance C g . 

We focus on the strong Coulomb blockade regime, assuming that the number of electrons which can be added to the 
dot is restricted to one. The physical features revealed in the corresponding bias window can be easily extended to a 
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FIG. 1: Schematic drawing of a carbon nanotube suspended over a trench. A bias voltage is applied between the STM tip 
(source) and substrate (drain). The inflated portion in the suspended portion of the nanotube illustrates the radial breathing 
mode, h is the tip-nanotube separation at rest. 



full span of many Coulomb blockade peaks, as well a several orbital levels as in Refs. iTlSl The extra-charge electron 
state is locally coupled to a phonon mode: in Refs. 00 this mode is identified with the radial breathing mode (RBM) 
of the carbon nanotub o?! 1 ? and, in Ref. 0, with the longitudinal stretching mode. In view of the phonon-assisted 
tunneling processes, the above-mentioned Coulomb blockade regime implies that the charging energy of the dot is 
assumed to be infinitely large compared to the relevant energy scale determined (at low temperature) by the phonon 
energy f2 and the bias voltage. In the absence of phonons, electron transport would occur only through a single 
electron state on the dot. 

In the Coulomb blockade regime, far from the Kondo regime, spin degrees of freedom are neglected. The Hamiltonian 
of the system is written as follows: 

H = Hq + Hleads + V , (1) 

where 

H = (e- gx)dU+Qtfb , x^b+tf, (2) 



H 



leads 



E 

jk 



(3) 



V = Y i T j (x)c] k d + h.c. . (4) 

jk 

Here the operator d (d)) annihilates (creates) an electron on a single dot level of energy e; similarly, Cjk (ct fc ) annihilates 
(creates) an electron with momentum k and energy £jk in the jth lead. The RBM is linearly coupled (with the coupling 
energy g) to the electric charge on the dot; the RBM excitations are annihilated (created) by b Qv). In Eq. Tj is 
the energy associated with the tunneling coupling to the dot; for simplicity, Tj is assumed to be energy-independent 
(constant density of states in the metallic leads), but we take into account the dependence of the tip-tube tunneling 
matrix amplitude on the boson coordinate x: due to the "breathing" motion of the tube, the tip-tube tunneling 
distance deviates from its equilibrium value. Explicitly, we assume an exponential ^-dependence: 

T l (x)=T lo e- sx , (5) 

where s is determined by the ratio of the amplitude of the zero-point RBM oscillations to the electronic tunneling 
length (sa 0.5A) which characterizes the tunnel barrier between the STM tip and the nanotube. Such position 
dependent amplitudes have been introduced in the context of nano-mechanical electronic devices^Li^SiM, where one 
refers to the "shuttle" mechanism as the central region oscillates between the two electrodes. On the other hand, in 
our situation the tube-substrate tunneling matrix amplitude is a;- independent: T r (x) = Tq. For s ^ 0, we refer here 
to this position dependent tunneling Hamiltonian as the "half-shuttle" mechanism, because only one of the tunneling 
amplitude (left) is modified by the position. 

In the model Hamiltonian (Q, formally describing the phonon-assisted tunneling of otherwise noninteracting elec- 
trons, the charging effects are taken into account via the bias-voltage dependence of the position of the dot level (e) 



4 



with respect to the chemical potentials of the leads, fii >r . A gate voltage V g and gate capacitance C g can be included. 
In our analysis, e includes the change in the charging energy of the dot when one extra electron is added to the dot, 
and which is obtained from the electrostatic energy consideration: 

e - e = 2E c (n x + 1/2) + e</> , (6) 

where eo is the "bare" energy level of the electron level in the dot, Ec — e 2 /2C is the charging energy and en x is a 
background (fractional) charge of the dot. Furthermore, the potential of the dot, <fr, is given by 

<t>= r ' • r ^ :: > ( 7 ) 



Cj = Cj/C , Cg = Cg/C , C = Cl + C r + Cg . (8) 

In typical experiments C g Cj . It is worth noticing that in the "floating-level" geometry which is considered in part 
of this work (e. g. without the gate electrode, C g = 0), the fractional charge is not fixed but it can instead be affected 
by changes in the capacitances of the junctions^ (that is, n x may depend on the tip-tube separation). 

Assuming weak coupling to the leads, it is convenient to eliminate the electron-phonon coupling in TIq by the 
unitary transformation H = A'7iA, with 

A = e -' iapdid , p = -i(b-tf) , a = g/Q, (9) 

In the rotated basis, the electron state in the dot becomes "dressed" with phonons, forming a small polaron. This 
results in the "polaron shift" of the dot level energy, e = e — g 2 / 17, and a renormalization of the dot-lead tunneling 
coupling. In the polaron representation, the Hamiltonian reads: 

H = 7Y + V + Hleads , (10) 

Ho = erfd + ntfb (11) 

V = Y,T j {x)e- ia »c\ k d+h.c.. (12) 

jk 

In Eq. H12J) . we have used the fact that 

tfTi{x)dk = Ti(x + 2adU)de~ iap = %{x) d . (13) 



III. RATE EQUATIONS 

Based on the assumption that the leads are in thermal equilibrium at given chemical potentials (fii and /i r ), 
independently of the state of the dot, one can derive a kinetic equation for the reduced density matrix by tracing 
out the electrode degrees of freedom from the total density matrix. Such equations have been used by several 
authorS ) 1:L i 13 i 16 i 32 but for sake of completeness we provide here a full derivation. 

In the polaron representation, the reduced density matrix of the phonon-coupled dot is given by 



R(t) 



-itH 



R v {t) 



Rv(t) = TrieadsPv(t) 



(14) 



where pv{t) is the total density matrix in the interaction picture (with respect to the tunneling coupling to the leads), 
which obeys the operator equation of motion: 



id t p v (t) = [V(t),p v (t) 
V(t) = e i *('^o+ W! <! c ' ds ) y e ^' t ( f 'o+'Hi e aj. 



(15) 
(16) 



Using the integral form of Eq. (|15|l , 



Pv{i) 



i I dr 
lo 



V(t),p v (t) +p v (0), 



(17) 
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substituting Eq. (|17fl into Eq. l|15fl , an d taking the trace over the electronic degrees of freedom of the leads with the 
above mentioned assumption of f>v(t) — Pleads ® ^vWi where p^ ads is the equilibrium density matrix of the leads, 
we arrive at the integro-differential equation of motion for the reduced density matrix in the polaron representation: 

d t Rv{t) = - 2rrv J J dTel " [ F /( T ) {^JW^K* ^ T ) dUR v (t - r) - %{t - r)dR v {t - r) St}{t)} 

3=1, r 

+ F<(t) { R v (t - t) dd)%{t - r)T/(i) - f]{t) S R v {t -T)dfj{t- r)}] + h.c. , (18) 
where Vj is the density of states in the jth electrode, 

f j {r)=T j [{x{r)]e- i ^ , (19) 

with 

x(t) = be~ inT +h.c. , (20) 
p( T ) = ~ibe~ mT + h.c. , (21) 

and 



5( T ) t P.v. 



/3sinh (ttt/jS) 



(22) 



(where "P.v." stands for the principal value). From the above expressions for the kernels i^ >:< (r), it follows that the 
relevant retardation time is of the order of the inverse temperature, (3. Assuming temperature to be high compared 
to the tunneling rates, /3 _1 3> I'KvfT^ and considering the long-time behavior of the reduced density matrix, t 3> j3, 
we construct an asymptotic solution of Eq. H18|) by using an improved perturbation expansion^: 

R v (t) = R(Xt) + A (t) + o(A) . (23) 

In Eq. H23|) , A is a formal perturbation parameter that reflects a weakness of tunnel coupling and must be formally 
attributed to the right-hand side of Eq. I|18|) : 

\d z R(z) + \d t R (1) {t) + o{\) = -\J2 27rv 3 J dTelir F>(t) {f}(t)f 3 (t-T)dUR{z-\T) + ... , (24) 

j=l,r 

with z — Xt. In the first-order equation in A, retardation effects are neglected (the Markovian approximation), and 
in the above equation, we can replace J Q dr ... — > J °° dr ... . By constraction, R is, on the time scale /3, a slowly 
evolving part of the reduced density matrix, while the rapidly oscillating (with frequencies determined by a multiple 
integer of fl) part of the right-hand side of Eq. (|24|) is absorbed into In what follows, we assume that the 

phonon frequency is high compared to the tunneling rates, f2 S> 2irvjT^ . This "anti-adiabatic" condition allows us to 
ncg lect i?W 

as being negligibly small compared to R. Under the above assumptions, Eq. I|24|) reduces to a differential 
operator equation for R{t). 

Although R is diagonal in the on-dot electron number (no coherence between subsequent tunneling events at long 
times), 

R(t) = (l-dU)R°(t) + dUR 1 (t) , (25) 

yet it is non-diagonal in the phonon number due to the presence of the displacement operators (e ±lap ) in the right- 
hand side of Eq. I|24|) . However, equations for the diagonal and off-diagonal parts of R(t) are decoupled from each 
other. As a result, in steady state (i — > +oo), which we are interested in, R becomes diagonal in both the electron 
and phonon numbers. 

Introducing the electron-phonon joint probabilities 

K(t) = (n\R*(t)\n) , (26) 

for i = {0, 1} additional electrons and n — {0, 1,2, ...} boson excitations on the dot, we obtain the system of rate 
equations for the joint probabilities, which in steady case can be written in the following form: 

m 

d t P l n - = -r> P r \ + J2 L< m P° , (27) 
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with the normalization condition, (i?^ + R\ 1=1. The charge/phonon transition rates are given by 

r ? i = ^ t L mn , r n = ^ t L mn , (28) 



L nm = E T i Ti,mn fj(^~^mn) , L> m = E T J lj,nm i 1 ~ M^M) • ( 2 9) 
j=l,r j=l,r 

In the above equations, fj(w) = [e^^ - ^^ + l] 1 is the Fermi distribution function of the jth electrode, fl nm = 
(n — to) fl, 

' • ^I'/li- (30) 

li, nm = (n\e- sx e- tap \m) , 7r ,„ m = (n| e -"* p |m) . (31) 
Explicitly, we have for the oscillator matrix elements: 



Min(n,m) m _^^^^ 



From the solution of Eq. I|27l) we can calculate the DC current and expectation values of the phonon observables. 
It should be noticed that in the original (non-polaron) representation, the resulting phonon density matrix is given 
by 

Rph = E ( P n I") ("I + Pn e- iap |n) (n| e^) . (33) 

n 

In general, R p h has contributions from both charge states i = 0, 1 and is non-diagonal in vibrational space (two shifted 
oscillators). 

In order to express the current in terms of the joint probabilities P„, we have to average the current operator taken 
in the polaron representation, 

Ij = ei J2 T o i x ) e_MP c )k d + h - c - » ( 34 ) 

k 

with the total density matrix 

p(t) = e -<n dot +H l<iads ) jt{H iot +H leaiB ) ; ( 35 ) 

where py (i) obeys Eq. I|17|l . Following the same steps and assumptions which we have made in the derivation of the 
rate equations, we obtain the average DC current flowing from the jth lead to the dot as (we set e = 1) 

^=E (Ff,n^-^j,nPn) > (36) 
n 

where are the partial relaxation rates contributing to r^ j> , Eq. l|28f) . from the jth electrode: 

T f,n = E T J rtnmfj(Z - {l nm ) , (37) 
m 

r ln = E r ^ i» c 1 - - ■ ( 38 ) 

m 

By virtue of the rate equations, we have conservation of the current, // = — I r . 
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IV. PHONON-ASSISTED TRANSPORT FOR T t < T r 



In calculating /(^-characteristics from Eqs. (|18|l and l|36|) . we will focus on the case of highly asymmetric two- 
junction model with T; <C T r . This corresponds to the typical experimental situation with STM measurements where 
a contact to the STM tip plays the role of the high resistive tunneling junction. For instance, in STM measurements on 
suspended nanotubes, the typical ratio Ru P /R S ub ~ 10 3 — 10 5 can be huge depending, in particular, on the tunneling 
distance between the STM-tip and the nanotube, h ~ 4 A typically. In the present Section, no gate is present (floating 
level) thus C g = 0. 

At the same time, according to the data in Ref. |8| obtained from spectroscopy measurements on suspended nan- 
otubes, the ratio Cu P /C su b = Ci/C r can be smaller as well as larger than unity depending on an effective length of 
the portion of the nanotube that is on the substrate, and which can be relatively short due to local defects induced 
by the edges of the trench. Thus, there is no dominating asymmetry in the capacitive coupling of the dot to the both 
electrodes. As a result, in the voltage-biased system, the effective position of the polaron level with respect to the 
chemical potentials of the leads is strongly affected by the ratio Ci/C r which is not negligibly small like r;/r r , or can 
even be larger than unitySi. It results that the polaron level (assuming neither gate-induced nor intrinsic shift of the 
polaron level at zero voltage) is not trivially stuck to the chemical potential of the far less resistive junction electrode, 
[i r . 

Considering the system of stationary rate equations, Eq. (|27|l . and writing explicitly expressions for the Fermi 
factors entering Eq. we obtain [see also Eqs. I©, Q)]: 

fi(e) = f(E-c r V) , fr(i) = f(E + ciV) , V = W -Mr, (39) 
where f(u>) = (e 0u + and 

E = e a - g 2 /n + 2E c (n x + l/2) (40) 

determines the position of the polaron level with respect to the lead chemical potentials at V = 0, considered as 
the reference level of zero energy. From Eq. (|39l) . the role of charging effects in the single-electron resonant level 
problem can be viewed as follows : when changing the bias voltage, the chemical potentials move in opposite directions 
(depending on the sign of V) with different "velocities" determined by the capacitance ratios, c r and c\ for the left 
and right electrode, respectively. The above picture is given in the "reference frame" of the polaron level, where its 
position is voltage- independent. For the floating-dot geometry, E may be strongly affected by a background charge, 
n x . In the presence of a gate electrode 3 ^, E can be controlled by the gate potential, and C in Eq. J7J), must also 
include the gate capacitance, with c\ + c r < 1. 

In this section, in order to make analytical progress on the prediction of the current transport, we formally assume 
temperature to be zero. Because the rate equations have been derived in the high-temperature approximation (T 3> 
r r ), the zero-temperature assumption is justified for the voltage-biased cases where the chemical potentials, /i;. r , are 
not very close to the phonon sidebands, so that thermally activated tunneling processes can be neglected. 

We first notice that in the case where the polaron level E lies beyond the bias- voltage window, the system of rate 
equations, Eq. i(27|) . reduces to the following: 

T< P° = T> P 1 (d}\ 

nm m ran n s V / 

which is satisfied for all phonon numbers m, n. The ground state solution of Eq. I|41|) is either _P* = <5„o<5,:i (i = 0, 1) 
for E < fj,^ r , or P^ — 5 n oSio for E > As a result, in this case the current is zero : no contribution to the current 
occurs from the phonon sidebands alone (cf. Ref. Il6|) . 

Eq. H41f) also works in the less trivial case when E lies within the bias-voltage window, but \^i,r — E\ < 51. For 
instance, for V = /i; — /i r > 0, the solution of Eq. I)41|l (valid for arbitrary Ti :T ) reads P„ — S n QP l , where 

P° = 1 - P 1 = T< L l° T> = , (42) 

tt = T ie - 2sa+s2 , (43) 

and the current is given by 

/ = P - a2 TlTr . (44) 

r* + r r 

For large a, the exponential prefactor in Eq. 144|) leads to a significant suppression of the current at low bias voltages 
(the so-called Franck-Condon blockadei^iSiii) . The role of s ("half-shuttle") will be discussed later. 
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At higher voltages, such that phonon sidebands enter the bias-voltage window, Eq. (|41() fails, and we have to 
go back to the full system of rate equations, Eq. i|27|). taking into account the total balance of probabilities for 
transitions between different states of the dot. This renders the problem of calculating phonon distribution functions, 
P^, and eventually the /(V)-characteristics, a rather complicated one to be solved analytically. However, for the case 
of high asymmetry in the tunneling couplings, it is possible to make some further analytical progress before passing 
to numerical calculations^. 

We now exploit the condition Ti <C r r for analyzing steady solutions of the rate equations. In the general case, 
when a bias voltage V is applied and the current is not identically zero (up to thermally activated processes), the 
polaron level is accompanied by N + M phonon sidebands captured in the bias-voltage window, where we introduce 
("int" stands for the integer part) 

N = wk(\m-E\/Q) , (45) 
M = mt(|/x r - E\/ti) . (46) 

Independently of the sign of V, the integer N (M) is defined as a number of phonon sidebands lying between the 
polaron level and the chemical potential of the electrode connected to the dot through the more (less) resistive junction. 
An example with N = 2 and M = 1 for V > is shown on Fig. [2J Due to high asymmetry T/ <C T r , only probabilities 
{P® n } with m £ [0,M] do not vanish in this limit. For the example on Fig. |2 the probabilities to have an electron 
on the dot with any n phonons, P^, are suppressed due to "fast" tunneling of the electron to the right electrode via 
n + 2 open channels; the probability to have the dot with an empty electron state but with the number of phonons 
A/ + 1 = 2 (and higher), P^ +1 , is also negligible due to the "fast" tunneling of an electron from the right electrode 
with absorbing two phonons. In other words, on a large time scale determined by T^ 1 , the polaron-hole states with 
and 1 phonon excitation are (quasi)steady states with respect to the tunneling coupling to the right electrode. Notice 
that in the asymmetric situation, the phonon number distributions are nearly the same in the polaron basis (the P„'s) 
and the original basis (R p h), as shown by Eq. I|33fl . 
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FIG. 2: Energy-level diagram for the case M = 1 and N = 2; the arrows show possible channels for electrons to tunnel 
onto/from the dot with changing (indicated by numbers) the phonon occupancy; fii — fi r = V 

The case V < can be treated in the same manner. We obtain that only probabilities {-P,^} with rn 6 [0, M] are 
not vanishing as r;/r r . As a result, with a good accuracy, controlled by the smallncss of r;/r r , the current flowing 
from the left to the right, I = Ii [see Eq. (ISBT) ]. can be written as 

M m+N 

I(V > 0) = P ° m Im , I,n ^^Y, limn , (47) 

m=0 n— 

wh ere £^ =0 P° + 0(r ; /r r ) = 1, and 

M m+N 

I(V<0)=-J2 P L Im, I m =r i Y J llnm , (48) 
m=0 n— 

with Era=0 p m + 0(r ( /r r ) = 1. The current is given by a sum of partial currents, I m , representing M conducting 
channels which arc distributed with the corresponding phonon occupation probabilities. The magnitude of each partial 
current depends on the previously defined number of phonon sidebands N captured between E and /i;. 

When varying V, the number of captured phonon sidebands, N and M, can change. In the voltage ranges where 
N and M do not change, we obtain plateaus in the /(V) characteristics. Assuming for concreteness V > 0, we 
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see from Eq. (|47|l that if V increases in such a way that only N changes (N can only increase with increasing V), 
the magnitude of each partial current, and hence the net current, will increase. Thus, we obtain here the positive 
differential conductance (PDC) behavior which is commonly computed and observed. On the contrary, if when 
increasing V, only M changes (i.e. increases), while N remains constant, then a new phonon-assisted channel is 
added. This leads to redistribution of phonon occupation probabilities between all open channels. Paradoxically, the 
net current may decrease, leading to a NDC behavior on 7(y)Il . 

For V > 0, assuming that due to the voltage increase V — > V + AV the number of phonon sidebands in the 
bias-voltage window has changed as {N = 0, M = 0} — > { N = 0, M = 1}, using Eq. (|47|) the variation of the current 
AI = I(V + AV) - I(V) becomes: 

AI = P °/ + P?h -h = -Pi tie- 01 " [ 1 - (a - sf - (l - a 2 + s 2 ) 2 ] . (49) 

Thus, for this particular situation, assuming for simplicity s = 0, the condition for the NDC is a < 1: large values of 
the electron-phonon coupling do not favor NDC. The magnitude of the negative step on /(V) increases with increasing 
phonon occupancy Pf. 

For V < 0, considering the same transition, {N = 0, M = 0} — > {N = 0, M = 1}, one finds: 

AI = -Pi f t e~ a2 \ 1 - (a + s) 2 - (1 - a 2 + s 2 ) 2 ] , (50) 

with the same NDC condition, a < 1, for s = 0. From the above expressions, we also notice that for positive 
(negative) bias voltage and non-zero but small half-shuttle s -C 1, the NDC effect increases (decreases) for the first 
phonon-assisted step of the /(V") characteristic. 

As far as the bias voltage increases, the current eventually saturates and does not change in practice. For large 
bias voltages V — > ±oo, the corresponding saturation currents 1^ are given by (up to corrections 0(r 2 /r r )): 

J« = r^P°(m|e- 2 -|m) , (51) 
lL } - -r, e- 4sQ £ PL (m|e- 2 Hm) . (52) 

m=0 

In the absence of the half-shuttle mechanism, we obtain that the saturation current does not depend on the phonon 
distribution, 7^ t = T;. For s ^ 0, considering equilibrated phonons, i.e., forcing — > 5 m o in the limit of strong 
phonon relaxation, we obtain the saturation currents: 

I^=r ie ^, lti=-e-*°"I%. (53) 

For non-equilibrated phonons, Eq. I|53|) is expected to be accurate in the limit q < Cr. Nevertheless, for relatively 
small a < 1, the asymmetry in the current saturation values due to the half-shuttle can be noticeable at low bias 
voltages for q ~ c r because of the fast saturation of the partial currents I m (fast convergency of the series in Eqs. 
(|47|> and gEJ}). 



V. NUMERICAL RESULTS AND DISCUSSION 



We now turn to the numerical solution of the rate equations, and explore the parameter space in order to find the 
signatures of NDC. After writting the system of the rate equations, Eq. (|27|l . in the matrix form, AP = 0, with 
P T = (P° , Pg 1 , P-f, Pi, ...J, the problem reduces to finding the zero eigenvalue of the matrix Aii& The typical maximum 
number of phonon states has been taken around 40, so that results of the calculation would not depend on this choice 
for the considered range of parameters. From numerically found solutions P, the current has been calculated in the 
left and right electrodes using Eq. I|36() . The achieved tolerance for = — I r was 10 -6 . 

Contrary to the previous analytical arguments, calculations are performed at finite temperature. For calculation 
purposes we assume a large asymmetry in the tunneling rates, r;/r r = 10~ 4 . Yet, most of our results also hold for 
moderate asymmetries. The current I is plotted in units of r;r r /(ri + T r ), the bias voltage V, the polaron level E 
and temperature T are given in units of f2. In subsection A, one considers the floating level case. We consider the 
effect of a given shift of the dot level (due to background or gate charges) in subsection B. In these two subsections 
we will disregard the half-shuttle (s = 0), which will be considered in the last subsection. 
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A. Non-shifted polaron level (E — 0) 

We start by considering the case where the capacitances surrounding the dot fully specify the position of the polaron 
level. At zero bias voltage, the polaron level is aligned with the chemical potentials of the electrodes, E = 0. In 
this case, as it immediately follows from the system of rate equations, Eq. I|27|) . we have the symmetry relation 
P®(V) = Pn(— V) and consequently I(—V) — —I(V). Notice that for the "floating- level" geometry the case of E = 
is an exceptional one, and must be viewed as a reference point. 

Fig. shows the I(V) characteristics for different a's assuming small c; = 0.1, so that the polaron level lies 
closer to [i r . In the experiment of Refs. I^ l30l this means that the nanotube has a large overlap with the conducting 
substrate. As long as M — 0, that is if V < Q/ci (this condition is satisfied on the bias- voltage range which is plotted), 
P® = S n o — 0(r//r, ), which results in the PDC behavior of I(V). The PDC steps correspond to N increasing by 
one each time V passes through a multiple integer of Q/c r . Explicitly, the current-step amplitude at V — n£l/c r 
corresponding to the current increase is given by 

rv 2n 

AJ< n > = T t e~ a . (54) 
n\ 

For smaller a's the current-step amplitude decreases faster with increasing V, which leads to the saturation of the 
current at lower voltages. 

Note that the first step (low bias voltage) is rounded. This feature is specific to the fact that the capacitances are 
asymmetric. Indeed, increasing the bias voltage from zero, the chemical potentials move away from the polaron level, 
but the chemical potential whose lead has the largest capacitance remains close to this level. In this event, thermally 
activated tunneling processes can be effective. When the bias voltage is further increased the phonon sidebands 
which contribute to new steps are those above the polaron level, and thus their corresponding steps are not thermally 
rounded. In the case where the asymmetry in the capacitances is opposite, which is shown in the next figure, the 
same reasoning applies. 




FIG. 3: The case q = 0.1 for a = 0.8 (solid), 1.5 (dashed) and 1.9 (dot-dashed). T = 0.01. 

FigJH shows the I(V) characteristics for different a's assuming a relatively large q = 0.9. In the experiment of 
Refs. Iffl this means that the nanotube has a small length, thus a small overlap with the conducting substrate, and 
thus a small capacitance compared to that of the tip-nanotube contact. Alternatively, a large ci could be achieved if 
the nanotube contains impurities or bends. 

In the bias- voltage window, the polaron level is closer to \ii than to fj, r . Within the bias- voltage range which is 
plotted, the current changes by steps when V passes through a multiple integer of fi/cj, as before. At these points, 
M increases by one; correspondingly, one more phonon-assisted channel is added. For a = 1.5 one still observes the 
PDC behavior, although the height of the current steps is strongly suppressed compared to the previous case on Fig. 
|3] For a < 1 NDC occurs in the first step, and with decreasing a, more NDC steps appear. For all curves, Eq. 
(14911 gives the current jump between the second and first plateaus on Fig. We have also performed calculations 
of the current voltage characteristics in this NDC regime while increasing gradually the temperature (not shown). 
Temperature tends to smooth the steps and to suppress NDC altogether when it becomes comparable to the spacing 
between the steps. 

At the same time that we monitor the electronic current, it is instructive to examine the phonon occupation numbers 
in order to quantify PDC or NDC behavior. The phonon distributions P^, which play the role of probabilities of 
open channels, at different bias voltages are shown on Fig. This figure demonstrates the increase of the number 
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^0.4 




FIG. 4: The case c; = 0.9 for a = 0.5 (solid), 0.9 (dotted) and 1.5 (dashed). T = 0.01. 

of phonon excitations with increasing M (See Eq. (|46|) ). It shows that the phonon number which are occupied is 
restricted by M — 1. Thus, the phonon distribution allows to understand the height of the PDC or the NDC steps. 
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FIG. 5: Phonon distribution, P„, for the case of Fig. 0]at the bias voltage V — 2 (circles) and V = 4 (triangles). 

In the two previous cases, two limits were considered : either c; <C c r or c; 3> c r . On the way to the more general 
situation when c/ and cy are comparable, let us now consider in more detail the capacitively symmetric case when 
q = c r = 0.5. The I(V) characteristics for a = 0.8 and 1.5 are shown on Fig. [S] For a = 1.5, the current shows a 
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FIG. 6: The case c; = 0.5 for a = 0.8 (solid) and 1.5 (dashed). T = 0.01. 

stepwise increased with essentially no NDC behavior. Note that because the capacitances are equilibrated, the first 
step does not display substantial rounding due to thermal effects. For a — 0.8, one observes a stepwise increase of the 
current, but the onset of each step is associated with a small spike, thus exhibiting NDC behavior. This NDC behavior 
is thus qualitatively different from the one previously observed for the case of asymmetric capacitances. The spikes 
bringing the NDC features to the I(V) characteristics are robust with respect to temperature. In Fig. [7|one sees that 
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these spikes are broadened with temperature, but NDC persists as long as the temperature is not comparable to the 
phonon frequency. On the I(V) for a = 1.5 (Fig. [BJ, the NDC-singularities appear only for higher voltages, however, 
the first phonon-assisted current step (at V — 2fi) is also deformed. 
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FIG. 7: Temperature dependence of the NDC-singularity for a = 0.8: T = 0.005 (solid), 0.01 (dashed) and 0.04 (dotted). 



At low temperature, it is interesting to make a zoom of a given step in order to see how the spikes evolve when the 
electron phonon coupling is varied. Fig. [SJshows these "singularities" on the second step of the /(V) characteristic for 
a wide range of a's (cf. Fig. |HJ). Note that the NDC behavior is more pronounced with increasing a = {0.5, 0.8, 1.1}. 
Indeed, in the limit of vanishing electron phonon coupling, one gets only one step due to the polaron (electron) level, 
and all NDC features are absent. There is therefore an optimal value of a which displays maximal NDC behavior. 
Further increasing a, NDC converts to the PDC case (the curves for a = {1.5, 1-9}). 
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FIG. 8: From up to down a = {0.5,0.8,1.1,1.5,1.9} for c ; = 0.5. T = 0.01. The dotted vertical lines around V = 2 
(corresponding to V± = 2 ± 0.05) show the width of the Fermi distribution function at given T. 



In this case of equal capacitances (q = cv), the polaron level lies in the middle of the bias- voltage window at all V, 
so that both N and M increase by one when V passes through a multiple integer of 20. In particular, at V = 2f2 the 
transition {N = M = 0} — > {N = M = 1} occurs : below V = 2fi the polaron level is the only one which lies within 
the bias voltage window, while immediately above V = 20, two phonon sidebands are simultaneously captured by 
the bias window. As a result, in the vicinity of V = 2f2 (how close depends precisely on temperature), we have a 
competition between the PDC- and "possible NDC"-type contributions to the current, corresponding to the processes 
of electron tunneling from the left electrode into the dot with emission and absorption of one phonon, respectively. 

In the transition region, V « 20, the current is given by (here we take into account the Fermi factors) 

1 

I(V)=J2 P ™( V ) T ™, (55) 

m=0 



m+1 

Ira = Tl lLnnf(-C r V - Q mn ) , 

71=0 



(56) 
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and P ° + Pf + 0(Ti/T r ) = 1. Explicitly, we have 



1 + a 2 f(fl - V/2) + P°(V) a 2 (-1 + a 2 + 1 [(1 - a 2 )(2 - a 2 ) - a 2 ] f(Sl - V/2)) 



(57) 



^From the numerics, it follows that up to V = V+ = 2fi + W > 0, where W is the half-width of the Fermi distribution 
function (see caption of Fig. [HJ), we have -PP(V) = 0, thus the current successively increases up to: 

I(V+) =T ie - a2 ( l + a 2 ) . (58) 

This effect is due to the high asymmetry in the tunneling rates T;^: the probability to have an empty dot with one 
phonon, Pf, remains negligible, 0(r;/r r ), until there is some probability to have thermally activated electrons of 
energy E — f2 coming from the right electrode. Only then, for V > V+, where /(fi — V/2) has practically reached 
unity, P®(V) starts deviating from zero. Thus, we have 

I(V) - I(V+) = P°(V)T ie - a2 'g(a) , V > V+ , (59) 



g(a) = a 2 (-1 + a 2 + \ [(1 - a 2 )(2 - a 2 ) - a 2 ]) 



(60) 



For V > V+, the V-dependence of the current is given by P®(V). With increasing V, P®{V) also increases and 




FIG. 9: g(a) function, Eq. CT . 



saturates to some quantity (depending on a) which is comparable with P" (the right-hand side plateau on Fig. [SJ. 
From numerics we find that at V = 2.5J1, for instance, for a = {0.5, 1.1, 1.9}, we correspondingly have P° w 1 — Pg = 
{0.35, 0.48, 0.38}. From Eq. Q59J) and Fig. Elwe see that the NDC plateau disappears for small a < 0.5; with increasing 
a, the NDC contribution from an open channel increases-decreases having the maximum at a ~ 1.15, and for a > 1.4 
the NDC contribution converts to the PDC one (cf. Fig. [BJ. 

Up to now, we have considered only strong asymmetric capacitances and equal capacitances. A common feature of 
these choices is that the current voltage characteristics displays steps (with or without NDC behavior) whose spacing 
in voltage is essentially always the same at relatively low voltages. Indeed, for the asymmetric capacitances plots 
each step is separated by (approximately) f2, while for q = c r , this spacing is 20. If one was to plot the differential 
conductance, as in the experiment of Ref. [jj, one would obtain a periodic series of peaks. What happens to this 
periodicity in an intermediate situation where the capacitances are comparable, but not equal ? Fig. 1101 shows 
current- voltage characteristics with two steps of different length for the case c r = 1.5c;. The number of phonon 
sidebands above the polaron level, N, increases by one with the voltage period fl/c r = (5/3)0; for these steps we 
always have a PDC behavior as it is demonstrated in the previous section. At the same time, the number of phonon 
sidebands below the polaron level M — > M + 1 with the voltage period O/c; = 2.5 57 : one more channel becomes 
open which yields the possibility of NDC. The larger a, the more PDC-steps at low bias voltages are observed. 



B. Shifted polaron level (E / 0) 

In this subsection, we consider the more general situation where, due to gate or background charges, the polaron 
level E is shifted by a constant value with respect to the chemical potentials of the leads at zero bias voltage. For 
E =/= 0, the symmetry I(V) = — 1(— V) is violated, moreover I(V) is not shifted in a trivial manner, meaning that 
it does not simply follow E. Fig. ^2 shows two current- voltage characteristics for different values of E. As it was 
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FIG. 10: a = 1.9 (solid), 1.5 (dashed) and 0.8 (dot-dashed) for c ; = 0.4, E = 0, T = 0.01, s = 0. 



mentioned in Sec. IV, the current is zero when the polaron level lies beyond the bias voltage window. For negative 
V, the polaron level becomes captured in this window at V = —E/ci — — 4E. At this point, the number of phonon 
subbands below the bare level - JV in this case - suddenly changes from N — to N = mt(c r E/ci$l) = int(3£'/0). In 
the curves of Fig. ^2 corresponding to E = 0.451 and 0.751, we observe a transition to JV = 1 and N — 2, respectively. 
This results in a fast saturation of the curves at negative voltages (the saturation is faster for the curve with larger 
E). 

Previously, we found that small ci leads to PDC for the first phonon assisted peaks. Here, although c; = 0.25 is 
relatively small in Fig. El f° r positive V, we observe NDC behavior already for the first phonon-assisted step on I(V) 
(around V = 251) for E = 0.7. This can be explained as follows. At small positive voltages, we first have a zero-current 
plateau until the polaron level becomes captured in the bias- voltage window (JV = M = 0) at Vq = E/c r . This leads to 
the first PDC step on the I(V) curve at Vq = (4/3).E. Then the second plateau persists until V\ = (51 — E)/ci < 251, 
which corresponds to the phonon sideband level E — 51 becoming captured in the bias window (N = 0, M = 1). 
Because a = 0.8 is less than the critical value 1 (see Sec. IV), this transition gives the NDC step. Next a plateau 
persists until, V2 — (£1 + E)/c r when the phonon sideband level E + £l is captured in the bias window (JV = 1, M = 1) 
resulting in a PDC step on I(V). The conditions to observe NDC for the first phonon-assisted step are therefore: 

c r Sl > E> (c r - cj)n > , a < 1 . (61) 

The width of the NDC plateau is V 2 -Vi = (E- (c r - /c 4 c r . Eq. J5TJ} is fulfilled for the case E = 0.751 on Fig. 
1111 but it fails for E = 0.451, where we have the PDC behavior of the first phonon-assisted step on I(V), for V > 0. 
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FIG. 11: E = 0.4 (dashed) and E = 0.7 (solid) for a = 0.8, T = 0.01, cj = 0.25. 

Fig. 1121 shows the temperature dependence of the NDC step for E = 0.751. With increasing temperature T, the 
phonon-assisted channel opens for higher voltages when the one-phonon excited state of the empty dot becomes 
unreachable for thermally activated electrons from the right electrode at energy E — 51 + W, (as before, W is the 
half-width of the Fermi function). As a result, the NDC region with a decreased slope gradually shifts towards the 
next step on I(V). 

Fig. 1 131 shows the opposite situation: the NDC appears this time with increasing temperature. At low temperature, 
when the Fermi half-width W is negligibly small, we only observe PDC behavior for the first step on I(V) at Vq — E/c r , 
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FIG. 12: The case a = 0.8, q = 0.25, E = 0.7 at different temperatures: T = 0.01 (thin solid), 0.02 (dashed) and 0.04 (bold 
solid) . 



where we suddenly have a transition to {N — 0,M = 1}. At this threshold, both the polaron level and the first 
phonon subband below it belong to the bias voltage window. Then we obtain a plateau in the /(V) curve up to 
V\ = (E + fi)/c r , where the next PDC step is associated with the increase of N by one. However, if the temperature 
is increased so that W becomes comparable with Q — E, now at V — Vq the occupation probability Pf will be 
negligibly small because of thermally activated electrons tunneling from the right electrode. Pf starts developing only 
at V > (fi - E + W)/ci w 1.217 leading to NDC behavior. Here we have used that at T = 0.04fi, the half-width of 
the Fermi function is about 0.2S7. At much higher temperature everything is of course washed out (not shown). 




FIG. 13: The NDC induced by temperature for the case a = 0.8, c ; = 0.25, E = 0.9 : T = 0.005 (dashed) and 0.04 (solid). 

We conclude this subsection by some comments on the role of the asymmetry in the tunneling rates T's on current 
transport. Fig. 1141 shows I(V) characteristics for the different cases of junction asymmetry. As it is seen from the 
plots, the NDC steps in the cases of high and moderate asymmetry (r//r r — 0.01 and 0.1, respectively) turn into 
PDC steps in the symmetric case (Ti = T r ). We also notice that the NDC steps are slightly shifted with respect to 
their PDC counterparts in the symmetric case. This shift is associated with the finite width of the Fermi function, 
and it has been already been discussed in the two previous figures. The height of the phonon-assisted steps in the 




FIG. 14: For a = 0.8, T = 0.02, E = 0.25, c ; = 0.5: Fi/r r = 10" 2 (solid), 0.1 (dashed) and 1 (dotted). The inset shows the 
phonon distribution at the bias voltage V — 4: circles for F;/r r = 10 -2 and triangles for r;/r r = 1. 
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asymmetric case is noticeably large compared to the symmetric one. This is related to the phonon distribution, an 
example of which (at V — 4f2) is shown on Fig. 1141 for both cases (symmetric and asymmetric). In the asymmetric 
case, the number of phonons is restricted by M, while in the symmetric case the phonon distribution is more spread 
out. As a result, in the symmetric case, the phonon- assisted contribution to the current is weakened. 



C. Effect of the half-shuttle on I{V) 



In this subsection we discuss the role of the half-shuttle mechanism on I(V) characteristics. 

It is worth to be mentioned that in the absence of the half-shuttle (s = 0), the rate equations l|27|) are invariant 
under the transformation 



V 



-V , E 



■E, P" 



Pi 



(62) 



As a consequence, the current- voltage characteristics Ie(V) and I-e(V) corresponding to the cases of the polaron 
level shifted by E and — E, respectively, are related as follows: 



Ie(V) = -I-e(-V) . 
The presence of half-shuttle mechanism (s ^= 0) violates this symmetry. 

1 1 ' 



(63) 




FIG. 15: Asymmetry of I(V) due to "half-shuttle" only: s = 0.1 (solid) and s = (dashed). The case of a = 0.7, c ; = 0.25, 
E — 0, T = 0.04. The inset shows the differential conductance for s = 0.1. 




FIG. 16: Asymmetry of I(V) due to "half-shuttle": s = 0.1 (solid) and s = (dashed). The case of a = 0.7, c ; = 0.75, E = 0, 
T = 0.04. The inset shows the differential conductance for s = 0.1. 

Fig. 1151 shows I{V) for the case when q = 0.25 is relatively small and where the polaron level is not shifted, both 
of which are favorable for PDC. The current steps are suppressed at negative voltages (differential conductance peaks 
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in the inset of Fig. I15f) . while at positive bias voltage, the current steps have a tendency to increase. From Eqs. i|47[l . 
(|48|l we obtain that the current-step amplitude at V = nQ/ 1 c r corresponding to the current increase is given by 

A/ (n) = Tl e -**-2saW (a + sign(V) S ) 2n 

n\ 

In other words, the half-shuttle mechanism works in favor of the formation of the polaron state: the probability of 
phonon- assisted tunneling onto the dot from the left electrode is increased, while the phonon-assisted tunneling from 
the dot to the left electrode is decreased. Eq. Il5.'i[l for the saturation currents illustrates this fact. 

The case of relatively large ci = 0.75 is shown on Fig. El Like in the previous case, here we also observe suppression 
of the current steps at negative voltages. At positive voltages, the NDC steps become more pronounced compared to 
the PDC ones. Expressions for the one-phonon assisted steps are given by Eqs. (|49|l . (|50|l . 
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FIG. 17: The case a = 0.7, c; = 0.5, E = -0.8, T = 0.04, s = 0.1. 




FIG. 18: The case a = 0.7, c ; = 0.5, E = 0.8, T = 0.04, s = 0.1. 

We have provided concrete evidence that the half-shuttle mechanism alone, with reasonable values of the parameters, 
cannot produce NDC when c/ is small and when the bare polaron level is not shifted upwards. Nevertheless, based 
on the results of Figs. El an d El we conclude that in situation when NDC is present (shifted level or reversed 
capacitances), the addition of the half-shuttle tends to emphasize the NDC features. 



VI. CONCLUSION 



To summarize, we have provided an in-depth study of NDC behavior in molecular quantum dots or transistors. 
While this study covers a wide range of experimental parameters such as capacitances and electron-phonon coupling, 
here the focus was put on the case where the tunneling coupling from the dot to the source and drain leads is 
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asymmetrical. This choice was motivated by the experiment of Leroy and coworkersi^S where a nanotube is suspended 
over a trench, allowing nearly free vibrations of the phonon breathing mode. Electrons flow from an STM tip to the 
suspended region of the nanotube, and their tunneling amplitude is much smaller than that of the nanotube to the 
substrate. It turns out that this assumption allows us to quantify NDC behavior by analytical means. Yet, moderate 
asymmetries, often encountered with two poor metallic contacts, qualitatively display the same physics. 

We have used a microscopic approach modeling the molcculc/nanotube as a quantum dot, which is justified based 
on the early observation of Coulomb blockade behavior in carbon nanotubes. A single phonon mode couples on-site 
to the quantum dot, in two different way. First, as the principal mechanism, it couples to the electron density on the 
dot, as in a polaron model. Second, the half-shuttle coupling was introduced, with the motivation that it should be 
present in the experiment of Refs. l7H30l when the nanotube vibrates, the tunneling distance between the nanotube 
wall and the STM tip oscillates accordingly, but the amplitude of this motion is believed to be rather small (or the 
order of the zero point motion, which can be smaller than O.lA). 

With these ingredients, the polaron transformation eliminates the electron charge coupling to the phonon on the 
dot, and transfers this coupling to the tunneling Hamiltonian. At this point, we have chosen to describe the situation 
corresponding once again to the experiment of Rcf. I7l l3fl where the time scale associated with tunneling events between 
the dot and the leads is large compared to the temperature, so that electrons evacuated in the leads effectively loose 
their phase coherence. On the other hand, the quantum mechanical nature of electron-phonon dynamics within the 
dot (polaron formation) is fully captured. This happens to be the regime were most molecular electronics experiments 
are performed nowadays, as these experiments typically do not require dilution refrigerator technology. From the 
density matrix of the total system, a kinetic equation for the reduced density matrix (with the leads degrees of 
freedom integrated out) was derived, leading to master equations for the electron population and their associated 
phonon numbers. 

At temperature high enough for this rate equation is valid, but low enough that thermally activated tunneling 
from the phonon sideband can be neglected, analytical predictions were made for the current-voltage characteristics, 
thanks to our assumption of a highly asymmetrical molecule-lead configuration: the positive or negative differential 
conductance behavior depends both on the location of the bare polaron level and on the occupation of the associated 
phonon sidebands - above and below this level - which are included in the bias voltage window. 

Turning to numerics, attention was first drawn to the role of the relative capacitances to the dot. If the 
molecule/nanotube is in better contact with the lead which evacuates the electrons (the substrate), than with the 
one who injects electrons, it is plausible to believe that the capacitance of this junction will be larger than that 
of the injecting junction. In this case, no negative differential conductance is obtained. When increasing the bias 
voltage, the contribution of an increasing number of phonon sidebands leads to a standard staircase behavior in the 
current-voltage characteristics. Thermal effects tend to round off the first step. In the experiment of Ref. F^l30l 
the use of several nanotube samples with different lengths allowed to vary the relative capacitance of the injecting 
and of the tunneling contact. This provided a motivation to study the so-called "reversed capacitance" case, which 
according to our analytical prediction could justify the presence of an negative differential behavior. This was indeed 
observed numerically, but it was emphasized that intermediate phonon couplings are needed to observe NDC: "small" 
electron-phonon coupling leads to barely noticeable NDC, while "large" electron-phonon coupling suppresses NDC 
altogether, thus the need to use non-perturbative techniques to study NDC. 

Do these results disqualify the possibility of NDC if contacts where the capacitances between left and right are 
comparable ? The answer is no, although NDC could be more difficult to detect in this case. While the global shape 
of the current voltage characteristics displays the staircase structure typically attributed to PDC, the onset of each 
step is followed by a small spike, thus displaying NDC behavior. The amplitude of this peak can be optimized with 
phonon coupling. For slightly asymmetric capacitances, the same effects are observed, but the staircase structure is 
modified as it involves two distinct periodicities. 

Because background charges may shift the bare polaron level, and because of the possibility of using a voltage gate 
on the molecule/ dot we also studied the possibility of NDC in this case. This suggests that when the level is shifted 
upwards, NDC is possible even when the capacitance of the injecting junction is small. 

Finally, we have proposed that a half-shuttle mechanism may play a role in a STM experiment. It can be detected 
by the asymmetry of the current-voltage curves. Yet, in our opinion, this mechanism cannot account for NDC alone 
and only tends to increase the heights of PDC and NDC steps in the current-voltage characteristics. Notice that 
approaching the tip closer to the molecule should strongly enhance the half-shuttle mechanism. In our numerical 
calculations, we set s = 0.1 assuming the amplitude of the zero-point RBM oscillations to be of the order of 1pm. 
Estimating the zero-point amplitude, we have considered a single-wall carbon nanotube as an elastic hollow cylinder 
with diameter do ~ 1 nm, with wall thickness ~ 1 A, and with a length corresponding to the relevant breathing part 
~ 10 nm. Accordingly, the undeformed volume of the cylinder is Vq ~ lnm 3 . The clastic properties of the nanotube 
cylinder are characterized with the average Young's modulus Y ~ 1 TPa (see, for instance, Ref. l37f) . Assuming 
7zO as 10 meV (fi is the frequency of the RBM), one finds for the relative extension of the nanotube due to the 
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breathing, 

YV ) (65) 

Thus, the estimated diameter change, ~ Sx do, associated with the zero-point oscillations is of the order of 1 pm. 

Overall, we have provided a rather complete account on the possible occurrence of NDC behavior in molecular 
electronic transport, due to phonon-assisted transport. This phenomenon is due to a distribution of the total spectral 
weight of the injected electron on the polaron levels included in the bias window. Taking into account the correct 
microscopic expressions of the Franck-Condon factors associated to phonon transitions, we have demonstrated that 
NDC effects are a fingerprint of the polaronic nature of charge carriers. While we cannot rule out that NDC could 
exist in molecular system with symmetric molecular junction, it was not observed within this model. On the other 
hand, we have clearly shown that for asymmetric tunneling rates, there exist a wide range of parameters which lead 
to NDC. 

A striking result of this work is the fact that a simple counting of the NDC steps in the current voltage characteristic 
(or alternatively, a counting of the peaks in differential conductance) provides direct information about how many 
phonon subbands contribute to transport. The question remains about whether this NDC behavior can be observed 
in experiments. Our model could apply to Ref. |6j, if the tunnel couplings are indeed asymmetric. Concerning Ref. 
0, NDC occurs but only for the first phonon sideband. This can be due to several factors : first, as a general rule, 
the height of PDC or NDC steps tends to decrease with increasing phonon sidebands number. Furthermore, if the 
electron phonon coupling is weak, NDC is present but its manifestation is weak. Finite temperature tends to smear 
these steps, so in order to observe several peaks in NDC, lower temperatures would be required, keeping in mind that 
the tunneling contacts should not be too large in order to "avoid" any Kondo like regime. Another aspect which we 
have not considered in this work is the possibility of phonon damping. The population of the phonon subbands may 
rearrange due to the coupling with the environment of the molecul o 15 i 16 . The environment could include electron-hole 
pairs generated in the leads, which should be reduced because of our weak coupling assumption. On the other hand, 
for a "large" molecular system, the optical phonon mode (the breathing mode of the nanotube of Ref. |2J) is also 
accompanied by other phonon modes such as acoustic-like modes, which have typically lower energies, and which 
therefore can play the role of a "bath" . Because of the importance of the occupation of the phonon subbands for 
NDC in our present work, we expect that substantial coupling to this bath will tend to suppress NDC. 
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